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1 Introduction 

For a meaningful source reconstruction, an 
appropriate source model has to be chosen before 
solving the inverse problem. A Singular Value 
Decomposition (SVD) of the measured MEG / EEG 
data yields not only the number relevant of singular 
values but also the corresponding spatial patterns 
and time loadings. The main drawback of this 
approach is that these patterns and loadings are 
orthogonal to each other. Based on the field patterns 
obtained by the SVD analysis, an Independent 
Component Analysis (ICA) has been performed to 
reveal field patterns that are mutually independent. 
Using simulated MEG and EEG data, the results of 
the ICA method are compared to those obtained by 
performing SVD only. 

2 Methods 

2.1 SVD 

The observed data are stored in an [s X t ] matrix M 
that contains 5 channels with t time samples each. 
The Singular Value Decomposition of M yields [1]: 

m = irz v 1 (i) 

U is an [s X s] matrix containing s orthogonal field 
patterns stored column-wise, X is an [ s X s] diagonal 
matrix with the corresponding singular values 
(sorted in descending order), and V is of the 
dimension [t X s] holding the orthogonal time 
loadings. 

2.2 ICA 

A thorough description of the Independent 
Component Analysis can be found elsewhere [2,3]. 
In short terms, the ICA method interprets the 
measured data as a linear mixture of an s- 
dimensional random vector, whose time loadings, or 
“components”, are assumed to be mutually 
independent: 

M - AS (2) 

A is the so-called mixing matrix, and the matrix S 
contains the independent components (stored row¬ 
wise). In “basic” ICA, the number of independent 
components is equal to the number of observed 


channels, i.e. A and S are of the dimension [s X s] 
and [s X t\, respectively. 

Solving the ICA problem is to find a matrix W 
iteratively, so that the linear transformation of the 
observed data 

S = WM (3) 

yields components that are as mutually independent 
as possible. Suitable techniques for solving the ICA 
problem can be found e.g. in [4,5]. 

To speed up the iteration process, the observed data 
M can be decorrelated by a linear transformation 
(called “prewhitening” or “sphering”) [3]. We use 
SVD for decorrelation, since this also reveals the 
weights of the components. One can reduce the 
dimension of the problem by setting all singular 
values to zero that are below a certain limit. 
Assuming that the noise is estimated properly, all 
values below noise level are discarded. This allows 
us to apply the ICA only on the signal subspace, 
greatly decreasing computation time. Equation (3) is 
modified to the following: 

S = WV kA (4) 

Fred is a [c X t ] matrix of the sphered and reduced 
input data obtained by SVD, W and S are of the 
dimension [c Xc] and [c X /], respectively (c is the 
number of singular values above noise level). 

In contrast to the SVD, the components obtained by 
the ICA do not underlie a fixed order. To sort the 
independent components, we determine the column 
norm of the mixing matrix A (equation (2)). With 
equation (4), A can be expressed as 

A-LTL W~ x (5) 

When sorting, one has to permute both the columns 
of A and the rows of S. The column norms of A are 
the weights of the ICA, analogously to the SVD. 

2.3 Experimental setup 

To investigate the differences between SVD and 
ICA, an analytical head model, consisting of one 
sphere (MEG) or three concentric shells (EEG), is 
used. Four dipoles are placed inside the innermost 
sphere. 
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Figure 1: Position of the 148 magnetometers and the 
four simulated dipoles. 

The MEG setup contains 148 magnetometers as 
shown in Figure 1. For the EEG simulation, 81 
electrodes (extended 10/20 system, not shown) are 
used. Both data sets contain 256 time points, 
sampling rate is 1ms. Simulations are performed for 
strongly and weakly overlapping time loadings 
(Figure 2) with gaussian noise added to yield a 
signal-to-noise ratio of about 15. 

After noise estimation (e.g. by using a signal free 
time interval), the original data is transformed to 
SNR values. The dimensionality of the ICA problem 
is reduced by rejecting all (normalized) time courses 
with a corresponding singular value less than 0.9. 
Finally, the (non-orthogonal) spatial patterns 
obtained by ICA are extracted from the mixing 
matrix, and the components are sorted by their 
impact to the data. 

To evaluate the separation capabilities of both 
approaches quantitatively, both methods are used as 
a spatio-temporal filter. Each component found is 
isolated from the other three by either setting all but 
one columns of the mixing matrix (ICA), or all but 
one singular values (SVD) to zero. One fixed dipole 
is fitted to the remaining data (whole time range). 
The fit results are compared to the original locations 
and orientations of the simulated dipoles. Since the 
SVD yields orthogonal input data for the 
reconstruction, one expects no meaningful dipole 
positions or orientations. 

All computations were performed with the Curry 
V4.5 (Neuroscan) software package. 

3 Results 

Figures 3 through 6 show the results of the 
decomposition of the MEG data matrix. The first 
four spatial patterns and time loadings obtained by 
ICA are shown. 



Figure 2: Simulated data containing non- and 
strongly overlapping sources, respectively. 


Decomposing the EEG data matrix yields very 
similar results and is not shown here. As expected, 
the SVD is hardly capable of separating the 
simulated patterns or time loadings. Using this 
method for filtering leads to spurious reconstruction 
results for both MEG and EEG data, as shown in 
Figure 7 and 8. The location error increases for 
smaller singular values with less impact on the data. 
On the other hand, the ICA was able to separate the 
non-overlapping dipoles nearly perfectly (the 
location error is below ) and achieved quite good 
results in the case of strongly overlapping dipoles. 
Due to the reduction of the data space to four 
“channels”, the ICA algorithm took less than one 
second on a standard PC to find all four 
components. 

4 Discussion 

The ICA can be an additional useful tool in 
preprocessing MEG / EEG data for source analysis. 
It is important to keep the following topics in mind: 

1. To determine the dimension of the ICA problem 
correctly, special care must be taken to proper 
noise estimation. If the number of included 
components is too low, no adequate separation 
will result. Searching for too many components 
will at least slow down the iteration. Moreover, 
trying to separate “sources” that are not truly 
independent (e.g. additional noise) could make 
the iteration fail. 






Figure 3: Spatial patterns obtained from a) ICA b) 
SVD, four simulated dipoles with no time overlap 
(MEG data). 


Figure 5: Spatial patterns obtained from a) ICA b) 
SVD, four simulated dipoles with strong time 
overlap (MEG data) 


ICA Time Courses: 1 to 4 PCA (SVD) Time Courses: 1 to 4 ICA Time Courses: 1 to 4 PCA (SVD) Time Courses: 1 to 4 




Figure 4: Time loadings from the four simulated Figure 6: Time loadings from the four simulated 

dipoles with no time overlap a) ICA b) SVD. dipoles with strong time overlap a) ICA b) SVD. 


2. Since the ICA is based on statistical properties 
of the observed data, a reasonable large number 
of time samples (at least some ten points) must 
be taken into account. On the other hand, 
computation time increases linearly with the 
number of samples. 


Only in the case where one strong component 
dominates the data, the first pattern and time loading 
obtained by the SVD could represent a meaningful 
result. 
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Figure 7: Reconstruction results (in mm) of single- 
dipole-fits applied on the simulated MEG data and 
using ICA and SVD as a spatio-temporal filter. The 
locations and orientations of the true dipoles are 
nearly identical to the ( TCA (no overlap)” data 
group and are therefore omitted for better legibility. 
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Figure 8: Reconstruction results (in mm) of single- 
dipole-fits applied on the simulated EEG data and 
using ICA and SVD as a spatio-temporal filter. The 
locations and orientations of the true dipoles are 
nearly identical to the i( ICA (no overlap)” data 
group and are therefore omitted for better legibility. 
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